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We study the anisotropic two-dimensional Hubbard model at and near half filling within a func- 
tional renormalization group method, focusing on the structure of momentum-dependent correla- 
tions which grow strongly upon approaching a critical temperature from above. We find that a 
finite nearest-neighbor interchain hopping is not sufficient to introduce a substantial momentum 
dependence of single-particle properties along the Fermi surface. However, when a sufficiently large 
second-nearest neighbor inter-chain hopping is introduced, the system is frustrated and we observe 
the appearance of so-called "hot spots" , specific points on the Fermi surface around which scattering 
becomes particularly strong. We compare our results with other studies on quasi-one-dimensional 
systems. 
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I. INTRODUCTION 

The one-band Hubbard model serves as a mini- 
mal model for various correlated electron systems, 
since it is capable of capturing a number of non- 
trivial phenomena which are due to the interplay 
between kinetic and potential energy. In one di- 
mension, numerous theoretical methods are avail- 
able which have led to a thorough understanding 
of the low-energy physics*! In higher dimensions, 
however, rigorous statements are scarce and many 
controversies remain. It is therefore natural to 
ask the question how the cross-over from one to 
two or higher dimensions takes place. Further- 
more, since the discovery of quasi-one-dimensional 
organic conductors and superconductors we have 
access to materials which are realizations of this 
physical situation. In these compounds many in- 
teresting observations were made during the last 
two decades, some of which are still calling for a 
conclusive theoretical description^ 

In this work we will consider the evolution of a 
model system near half-filling, upon increasing the 
dimensionality via an increase in the perpendicular 
kinetic-energy coupling between one-dimensional 
chains. To tackle this question, we employ a 
functional renormalization group (fRG) technique, 
which provides a rigorous framework for the com- 
putation of low-energy properties starting from a 
microscopic modeli- While this method reduces to 
the well-known g-ology RG in one dimension, it has 
been applied successfully to the two-dimensional 
case, where the angular dependence of the coupling 
function along the Fermi surface needs to be taken 
into account^ By allowing for anisotropic hopping 
parameters we are in principle able to access the 
complete region from the one-dimensional to the 
two-dimensional case. In this work we will fix the 



degree of anisotropy and study the behavior of the 
system as a function of other parameters. 

The dimensional cross-over in quasi-one- 
dimcnsional systems (and the intimately related 
phenomenon of "deconfmement" , i.e the Mott 
insulator to metal transition induced by increas- 
ing the inter-chain kinetic energy) have been 
investigated in recent years within several theo- 
retical approaches. There is consensus on some 
aspects, but some disagreements between these 
studies do remain (mainly due to the different 
theoretical tools which have been employed, and 
the different regimes of parameters which have 
been investigated). The issue as such was raised 
by experimental studies on Beechgard salts. In 
particular, the optical spectroscopy experiments 
of Vescoli et al. revealed an insulator-to-metal 
transition as a function of increasing interchain 
hopping parameter, which changes with chemical 
composition 5 . At the same time, Bourbonnais and 
Jerome discussed these results within a scenario 
where a one-dimensional Mott insulator evolves 
into a metallic state when the inter-chain hopping 
reaches the order of the Mott gap£. 

Subsequently, several model calculations where 
made to substantiate and verify this concept. Bicr- 
mann et al. employed an extension of dynami- 
cal mean- field theory (DMFT), the so-called chain- 
DMFT, which replaces the original problem by a 
chain self-consistently coupled to a bath, while tak- 
ing into account the full intra-chain momentum de- 
pendence (0, see also ||). They do indeed find 
a transition from an insulating to a metallic state 
when the interchain hopping is increased, as well as 
a crossover from a Luttinger liquid to a Fermi liq- 
uid at fixed interchain hopping, when the temper- 
ature is decreased. Essler and Tsvelik considered 
this problem starting from a one-dimensional Mott 
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insulator and using a resummed expansion in the 
inter-chain hopping 9,10 . Using this approach, they 
suggested that the metallic phase does not develop 
immediately with a large Fermi surface resembling 
the non-interacting one. In contrast, close enough 
to the Mott insulating phase, Fermi surface pock- 
ets appear in specific locations, while large parts 
of the would-be Fermi surface remain gapped due 
to the influence of the one-dimensional Mott gap. 
Since they use a particular dispersion in the inter- 
chain direction which is particle-hole symmetric, 
they eventually find that this pocket Fermi liquid 
becomes unstable towards an ordered state, how- 
ever at much lower temperatures as compared to 
the Mott gap. The location of these pockets is such 
that the neighbourhood of the point k a — kb = 7r/2 
(with k a the momentum along the chain direction 
and kb perpendicular to it) is gapped out and not 
part of the Fermi surface. This point thus corre- 
sponds to a "hot spot", at which the scattering 
rate is very large (and can even lead to a complete 
suppression of quasiparticles at this point). While 
the chain-DMFT studies of Ref. 0] did not observe 
this phenomenon (perhaps because of the range of 
coupling or temperature), more recent studies of 
an anisotropic spinless model using chain-DMFT 
did observe a partial destruction of the Fermi sur- 
face with hot spots at the same location^. 

From the weak-coupling standpoint, early calcu- 
lations suggested the occurrence of hot spots from 
a simple perturbative calculations of the scatter- 
ing rate 11 . Within a renormalization group treat- 
ment, Duprat and Bourbonnais found that for a 
finite and fixed value of the interchain hopping 
the influence of strong spin-density-wave fluctua- 
tions can lead to anisotropic scattering rates along 
a quasi-one-dimensional Fermi surface, leading to 
the emergence of hot spotsi^. However, the loca- 
tions of cold and hot regions are exactly exchanged 
with respect to the findings by Essler and Tsvelik, 
with the hot spots found at kb = and kb = ±7r 
in Ref. ^ should be noted that these results 
were obtained for a system away from half filling, 
meaning that umklapp processes are suppressed in 
the RG treatment. In the present work, we fo- 
cus on the half-filled case or its immediate vicinity, 
and use a functional RG technique which does take 
umklapp processes into account. 

Thus, there exists a consensus that "hot spots" 
might form in quasi one-dimensional systems, but 
the various treatments do not agree on their lo- 
cation. It is tempting to speculate that this sim- 
ply reflects the different locations of these hot re- 
gions in the weak and strong coupling limits. At 
any rate, there are some compelling experimen- 
tal indications for strongly anisotropic scattering 
rates in quasi one-dimensional organic conductors, 
as pointed out early on in Ref. 14] and dicussed 



further in the conclusion. 



II. MODEL 

We consider the one-band Hubbard model: 

with a local repulsion U > and hopping ampli- 
tudes ijj' = — t a between nearest neighbors in the 
a-direction along the chains, tjy = —% between 
nearest neighbors in the b-direction perpendicu- 
lar to the chains, and tjy = —t b between second- 
nearest neighbors in the b-direction. The corre- 
sponding dispersion relation reads 

£k = —2t a cos k a — 2tb cos kb — 2t b cos 2kb- (2) 

At half-filling and t b = the non-interacting Fermi 
surface is perfectly nested. The introduction of a 
finite t b will allow us to study the effects of de- 
viation from this perfect nesting condition. It is 
important to note that we do not linearise the dis- 
persion in the chain direction, as is commonly done 
in RG treatments originating from the Id g-ology 
setup Therefore, the Fermi surface is perfectly 
nested only at half filling. Away from half filling 
perfect nesting is destroyed, even without a finite 
value of t b . In the following all energies are given 
in terms of t a which we set to unity. 

III. METHOD 

In direct analogy to the technique applied in ref- 
erence 0| we use the Wick-ordered version of the 
fRG to compute the one-loop flow of the interac- 
tion vertex and the two-loop flow of the self energy, 
as depicted in Fig. 1. The internal lines without 
slash in the Feynman diagrams correspond to the 
bare propagator 

D , (k) _ e<A-M , (3) 

ik - 4 k 

where £k = £k — A* an( l A > is the cutoff; the lines 
with slash correspond to d\D A , which is propor- 
tional to 5 (A - |£ k |). 

We parametrize the interaction vertex T by 
its static values on a reduced number of 
points/patches on the non- interacting Fermi sur- 
face, as illustrated in Fig. 2. It is thus parametrised 
by a momentum-dependent singlet (triplet) com- 
ponent T s n\(Jci, k2i A^), where the fc, constitute a 
discrete set of momenta on the Fermi surface. We 
stress that this does not correspond to treating a 



FIG. 1: Flow equations for the self energy E and the 
two-particle vertex T A , respectively; the internal lines 
without slash correspond to the bare propagator D A , 
the lines with a slash to its A-derivative 8aD a . 
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FIG. 2: Patching scheme for N = 28. 



finite system, since internal integrations are done 
in the thermodynamic limit. 

In the present work, we do not directly compute 
the flow of the self-energy, but rather infer from the 
properties of the two-loop diagram what the influ- 
ence of a strongly renormalized vertex function on 
the self energy will be. In the case of the two- 
dimensional Hubbard model this calculation has 
been done explicitly^ giving us confidence with 
respect to this reasoning. 
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FIG. 3: Fermi surface and nesting vector for tt — 0.1 
and t b = 0, that is the case of perfect nesting (All 
energies in units of t a = 1). The picture on the right 
shows the warping of the Fermi surface near k a = -j. 



IV. RESULTS 

A. t b = - perfect nesting 

At t h = and /i = the non-interacting Fermi 
surface is perfectly nested with nesting vector Q = 
(■7T, 7r) and defines the so-called Umklapp surface 
as illustrated in Fig. 3. For tb — the problem 
reduces to the one-dimensional half-filled Hubbard 
model. In this case the divergence of Umklapp cou- 
plings at low energies signals the onset of the Mott 
insulating phased When we analyze the one-loop 
flow of the interaction vertex for finite values of the 
interchain hopping tb, we find that the divergence 
of Umklapp processes persists at all finite values of 
tb, which we chose to be tb = 0.1i a throughout this 
work. Namely the one-dimensional Umklapp cou- 
plings connect to two-dimensional Umklapp pro- 
cesses of the type (kp,k' F ) — > (kp + Q, k' F — Q) with 
momentum transfer Q = (71", tt). The crucial point 
is that this divergence is nearly perfectly homoge- 
neous along the Fermi surface. Due to the feedback 
of the interaction onto the self energy via the two- 
loop diagram, this implies that there are no iso- 
lated hot spots at which the dominant scattering 
processes are dominant compared to other regions 
on the Fermi surface. Instead, the whole Fermi 
surface is "hot". Along with Umklapp processes 
the interaction develops divergences in the Cooper 
channel, also owing to the importance to scattering 
with wave vector Q = The behavior of the 

coupling function illustrates this very clearly. In 
the left plot of figure 0] we display the singlet com- 
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FIG. 4: Left: Singlet Vertex in the Cooper channel on 
the Fermi surface as a function of kt and k' b . Right: 
Singlet Vertex in the Twin-Umklapp channel on the 
Fermi surface as a function of k^ and k b . Here tb = 
0.1, t' b = 0, T = 0.011, fj, = 0, U = 2.45 (All energies 
in units of t a — 1). In this case the Fermi surface is 
essentially perfectly nested. 



ponent of the interaction function Fs(k, — k, — k') 
on the Fermi surface as a function of kb and k' b at 
the end of the flow, in direct analogy to the anal- 
ysis presented in reference |l2j . The bare interac- 
tion is chosen in such a way that the temperature 
is slightly above the pairing temperature, at which 
the flow of the vertex function diverges. The in- 
teraction is homogeneously peaked along the lines 
k' b — n — kb and k' b = —ir — kb- We name these lines 
"Peierls lines" since they correspond to scattering 
processes in which a momentum Q is exchanged 
between the two incoming particles and Q is the 
generalization of 2kF in one dimension at half fill- 
ing and perfect nesting. 

We thus see that all points on the Fermi sur- 
face are equally strongly affected by strong correla- 
tions appearing in the Cooper channel, which may 



eventually lead to a phase where all one-particle 
states along the Fermi surface exhibit a pseudo- 
gap. Since we look at a system at half filling it is 
essential to consider Umklapp processes, as men- 
tioned above. We therefore extend the analysis 
offered in reference |l2j and show in figure 01 in the 
right plot the interaction function Ys{k b ',kP , ,k^) 
on the Fermi surface in the so-called twin-Umklapp 
channel, meaning that both incoming momenta are 
identical, and k b {k b ) are momenta corresponding 
to right (left) movers in the standard terminology 
familiar from the one-dimensional case. We see 
that along the lines — k b — ir and k b = k b + n 
the coupling function in this channel is also homo- 
geneously peaked, confirming the conclusion drawn 
on the basis of the Cooper channel. We recall that 
in contrast to the one-dimensional case the system 
will eventually undergo a transition into an an- 
tiferromagnetic state at zero temperature. Here, 
however, we are concerned with finite-temperature 
precursor effects, which may legitimately be com- 
pared. 

These results are distinct from those obtained 
in other RG calculations^. There, perfect nesting 
is artificially introduced due to a linearization of 
the dispersion in the chain direction, and Umklapp 
processes are neglected. The divergent couplings 
are then found only in the Peierls section of the 
Cooper channel for sufficiently small t b , leading to 
isolated hot spots at kb = and kb = ±7T. 



B. t b 7^ - effects of frustration 

1. fi = - half filling 

A finite second-nearest neighbor hopping t b will 
destroy the nesting condition for all wave vec- 
tors on the Fermi surface, except for a few spe- 
cial points. We set the chemical potential to 
fj, = and use tb = 0.1 and t b = O.lif,. In 
this case the non-interacting system essentially re- 
mains half filled and the frustrated Fermi surface 
intersects the Umklapp surface at kb = ±7r/4 and 
kb = ±37r/4, as shown in Fig. [3J Then, at ar- 
bitrarily low energies Umklapp scattering of the 
type (kF,kp) — > (&f + Q,&f — Q) is possible if 
and only if kp is located at the intersection be- 
tween the non-nested Fermi surface and the Umk- 
lapp surface. The resulting RG flow of the inter- 
action vertex shows a dominant divergence of the 
couplings corresponding to exactly these processes. 
This can best be seen in figure where the cou- 
pling function is shown in analogy to figure In 
contrast to the case t' = the coupling function 
in both channels shows a strongly peaked behavior 
which is not homogeneous along the Peierls lines, 
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FIG. 5: Fermi surface for 4 = 0.1 and t b — 0.1t b (red, 
dotted line) and for t b = (black line). The points at 
which the two intersect allow for Umklapp scattering 
of the type (&f, Icf) — > (fcp + Q, &f — Q) at arbitarily 
low energies and therefore emerge as hot spots in the 
RG flow. 

but which is peaked at the points where the Fermi 
surface intersects the Umklapp surface. This is 
reminiscent of a scenario in which there exist so- 
called hot spots, that is special points at which the 
scattering rate is particularly large or a pseudogap 
may appear in the spectral function, in analogy to 
the case of a two-dimensional system^ Note that 
for the frustrated system the Peierls lines defined 
above in the plots of the vertex function correspond 
to scattering processes with wave vector Q only at 
hot spots. Elsewhere the momentum transfer on 
the Fermi surface is incommensurate. 



2. fi / - slightly doped system 

Upon changing the chemical potential the hot 
spots mentioned above move along the Fermi sur- 
face. When /z is increased, the two points in each 
quadrant eventually merge until the Fermi surface 
touches the Umklapp surface at (±7r/2, ±7t/2). In 
figure H we show plots for the vertex function in 
analogy to figures 01 and El Indeed, the vertex 
function along the Peierls lines in both Umklapp 
and Cooper channel exhibits a strong increase in 
the diagonal region, corroborating the identifica- 
tion of the hot spots as originating from the points 
where the Fermi surface intersects the Umklapp 
surface. For even larger values of the chemical 
potential these points do not exist anymore and 
thus Umklapp processes will not feed back into the 
self energy. Similarly, upon decreasing \jl the eight 




FIG. 6: Same as figure^ but here t = 0.1, t' = 0.14, 
T = 0.011, fi = 0, U = 3.08. In this case perfect 
nesting is destroyed and the system is frustrated. 

hot spots move towards the axis and eventually 
merge to form four hot spots located at kb = 
and kb = tt. Once more the structure of the vertex 
function reflects this, as can be seen in figure |3 
although for the parameters chosen here the vari- 
ation along the Fermi surface is somewhat weaker. 
For small enough values of /x the hot spots will 
again disappear. 



V. DISCUSSION AND CONCLUSION 

In summary, we have studied a quasi one- 
dimensional model of coupled chains at and near 
half-filling, using an fRG technique and focusing 
on the appearance of "hot" regions on the Fermi 
surface. In the presence of perfect nesting, we have 
found that the whole Fermi surface is hot. In con- 
trast, in the presence of frustration (t' b ^ 0), iso- 
lated hot spots appear. The mechanisms for the 
formation of these hot spots is that the effective 
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FIG. 7: Same as Figure © but here fx = 0.02, U = 
3.126. In this case the Fermi surface touches the Umk- 
lapp surface at (7r/2, tt/2). 



couplings (vertex functions) in the various chan- 
nels become large in an anisotropic manner. The 
location of the hot spots corresponds to the in- 
tersection of the Umklapp surface with the Fermi 
surface. In general, there are eight hot spots, i.e. 
two in each quadrant of the Brillouin zone. Their 
precise location depends on the doping level and 
on the ratio t' b /tb- These results are perfectly con- 
sistent with previous weak-coupling fRG studies 
of the frustrated two-dimensional casein. The lo- 
cation of these hot spots do not agree however 
with a previous weak-coupling RG study of a quasi 
one-dimensional model |I2j. As we have seen, a 
proper treatment of Umklapp processes, which are 
a key ingredient to the mechanism described in the 
present work, is essential. Because Ref.|l2| was 
motivated by the strongly metallic regime, it did 
not explicitly include these processes in the RG 
treatment, besides a mere renormalisation of the 
forward scattering amplitude. 

There are naturally also some limitation to our 




FIG. 8: Same as Figure OH but here fx = -0.0197, 
C/ = 3.75. In this case the Fermi surface intersects the 
Umklapp surface near kb = and kb = 7T. 



fRG calculation. First, it is valid only in the weak- 
coupling regime. Second, the argument that a 
strongly peaked interaction can create hot spots 
on the Fermi surface relies on the low-dimensional 
properties of the model, the reason being that with 
increasing dimension the feedback of the interac- 
tion onto the one-electron self-energy via the two- 
loop diagram weakens and is eventually washed 
out. The method does take transverse Umklapp 
processes into account, which is essential to the 
main mechanism and observed features. This can 
already be achieved by a much simpler RPA calcu- 
lation. However, the fRG not only provides an ex- 
act starting point relying on rigorous statements, it 
also modifies the RPA results. The critical scales 
are lower and in contrast to RPA the properties 
in Cooper and Umklapp channels are different at 
equal momentum transfer. 

A different location of the hot spots was also 
found by strong-coupling techniques, such as the 
resummation of the expansion in the inter-chain 
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hopping of Ref. $Httl and the recent chain-DMFT 
treatment of Ref. |L| ■ This is less surprising, and 
it is tempting to speculate that the hot spot loca- 
tion may be determined by the regions in momen- 
tum space where the effective couplings are big in 
the weak coupling limit, while it is associated with 
regions in which the inter-chain kinetic energy is 
small in the strong coupling limit. Future studies 
at intermediate coupling are needed in order to elu- 
cidate this point and provide a consistent picture 
of how the location of the hot spots evolve from 
weak to strong coupling. 

The possibility that electron-electron Umklapp 
scattering may account for the emergence of hot 
spots along a quasi- Id Fermi surface was first sug- 
gested by Chaikin in order to explain magic an- 
gles in the magnetoresistance data of Beechgaard 
saltsii Here we have shown on the basis of a 
microscopic model and a functional renormaliza- 
tion group approach how such a situation may 
arise. Angular dependent magnetoresistance oscil- 



lation experiments only provide indirect evidence 
for the formation of hot spots however, through 
a momentum dependence of the scattering rate 
on the Fermi surface which requires theoretical 
modelisation. Obviously, direct angular-resolved 
spectroscopy experiments (e.g photoemission), al- 
though very difficult to perform on quasi one- 
dimensional organic conductors, would be highly 
desirable in order to probe these effects experimen- 
tally. 
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